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Abstract 

We present a generalization of the hyperspherical harmonic formalism to study systems made of 

<^ . quarks and antiquarks of the same flavor. This generalization is based on the symmetrization of 

^:^ i the A'^— body wave function with respect to the symmetric group using the Barnea and Novoselsky 

5_^ I algorithm. The formalism is applied to study four-quark systems by means of a constituent quark 

Qh' model successful in the description of the two- and three-quark systems. The results are compared 

"^ ' to those obtained by means of variational approaches. Our analysis shows that four-quark systems 

cn . with exotic 0^ and non-exotic 2~^~^ quantum numbers may be bound independently of the mass 

of the quark. 2^ and 1"' states become attractive only for larger mass of the quarks. 
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I. INTRODUCTION 

The understanding of few-body systems relies in our capability to design methods for 
finding an exact or approximate solution of the iV— body problem. In two-, three-, and four- 
body problems it is possible to obtain mathematically correct and computationally tractable 
equations such as the Lippmann-Schwinger, Faddeev and Yakubovsky equations describing 
exactly, for any assumed interaction between the particles, the motion of few-body systems 
[l|. However, the exact solution of these equations requires sophisticated techniques whose 
difficulty increases when increasing the number of particles. 

There are a countless number of examples of quantum-mechanical few-body systems, 
from few-electron quantum dots in solid state physics to constituent quarks in subnuclear 
physics. The intricate feature of the few-body systems is that they develop individual 
characters depending on the number of constituent particles. The most important cause 
of these differences are the correlated motion and the Pauli principle. This individuality 
requires specific methods for the solution of the few-body Schrodinger equation, approximate 
solutions assuming restricted model spaces (as for example the mean field approximation) 
faihng to describe these systems. 

The solution of any few-particle system may be found in a simple and unified approach. 
A recent widely used method is the stochastic variational j^, a variational approach where 
the trial wave function is generated by a random search on an adequate basis. An alternative 
powerful tool is an expansion of the trial wave function in terms of hyperspherical harmonic 
(HH) functions. The idea is to generalize the simplicity of the spherical harmonic expansion 
for the angular functions of a single particle motion to a system of particles by introducing a 
global length p, called the hyperradius, and a set of angles, Q. For the HH expansion method 
to be practical, the evaluation of the potential energy matrix elements must be feasible. The 
main difficulty of this method is to construct HH functions of proper symmetry for a system 
of identical particles. This is a difficult problem that may be overcome by means of the 
HH formalism based on the symmetrization of the A^— body wave function with respect 
to the symmetric group using the Barnea and Novoselsky algorithm [3|. This method, 
applied in nuclear physics for A^ < 7 y, has only been applied to quark physics for A^ = 3 
J5|. Therefore, its generalization would be ideally suited for the study of the properties of 
multiquark systems. 

During the last few years there has been a renewed interest on the possible existence 
of multiquark states, specially four- (two quarks and two antiquarks) and five-quark (four 
quarks and one antiquark), in the low-energy hadron spectroscopy. Theoretically, the possi- 
ble existence of four-quark bound states was already suggested thirty years ago, both in the 
light-quark sector by Jaffe p and in the heavy-quark sector by Iwasaki [7|. Experimentally, 
there were several analysis suggesting the existence of non-qq states in the high energy part 
of the charmonium spectrum ^] . The recent series of discoveries of new meson resonances 
whose properties do not fit into the predictions of the naive quark model, has reopened the 
interest on the possible role played by non-gg configurations in the hadron spectra [9|, |lO| . 
Among them one could mention the X(3872) [ll|, the y(4260) [l^ or the new D and Dg 
resonances |l3|. Out of the several interpretations proposed for these states, those based on 
four-quark states are majority |l4]. There seems to be a consensus that a four-quark system 
containing two-light and two- heavy quarks, qqQQ, is stable against dissociation into two 
mesons, qQ, if the ratio of the mass of the heavy to the light quark is large enough |l5lll6l|. 
While the conclusions agree for the qqbb system, there are more discrepancies about the qqcc 



system QHH- 

The possible existence of four-quark states in the low-lying meson spectra has been used 
to advance in the understanding of the light-scalar mesons, J^^ = 0~^~^. In this case, four- 
quark states have been justified to coexist with qq states in the energy region below 2 GeV 
because they can couple to 0++ without orbital excitation p , and therefore conventional qq 
states are expected to mix with four-quark states to yield physical mesons [10|, |20| . 



The four-heavy quark states have not received as much attention in the last years as 
the ones containing light quarks. Iwasaki |3|, based on the constituent quark mass value 
obtained within a string model for hadrons, argued that four-charm quark states could exist 
in the 6 GeV energy region. Using a variational method with gaussian trial radial wave 
functions and two-body potentials based on the exchange of color octets between quarks, 
Ader et al. J^] obtained non stable QQQQ states. However, using a potential derived from 
the MIT bag model in the Born-Oppenheimer approximation the same authors concluded 
that the cccc state was bound by 35 MeV. Heller and Tjon [la| considered also the MIT 
bag model improving the Born-Oppenheimer approximation with a correct treatment of the 
kinetic energy, non stable QQQQ states were found. In Ref. ||21| L = QQQQ states were 
analyzed in the framework of a chromomagnetic model where only a constant hyperfine 
potential, not depending on radial coordinates, was retained. With these assumptions no 



bound QQQQ states were found. A similar conclusion was obtained in Ref. jl8[ using the 
Bhaduri potential and solving variationally the four-body problem in a harmonic oscillator 
basis. Recently, Lloyd and Vary [22] investigated the cccc system using a nonrelativistic 
hamiltonian inspired by the one-gluon exchange potential diagonalizing the hamiltonian in 
a harmonic oscillator basis, obtaining several close-lying bound states. 

The discussion above illustrates that the theoretical predictions for the existence of four- 
quark systems differ depending basically on the method used to solve the four-body problem 
and the interaction employed. It is our aim in this work to make a general study of four- 
quark systems of identical flavor in an exact way. For this purpose we will generalize the HH 
method J23|, widely used in traditional nuclear physics for the study of few-body nuclei, to 
study four-quark systems. There are two main difficulties, first the simultaneous treatment 
of particles and antiparticles, and second the additional color degree of freedom. As a test 
of our formalism we will recover some results present in the literature. Finally, we will make 
a general study of four-quark systems of identical flavor by means of the constituent quark 
model of Ref. 2J| that provides with a realistic framework describing in a correct way the 



general aspects of the meson and baryon spectra. 

The paper is organized as follows. In the next section the formalism necessary to build 
a color singlet wave function with well-defined parity and C— parity quantum numbers is 
described and discussed. In Sect. IIIII we introduce the general features of the constituent 
quark model used. In Sect. II VI we present and analyze the results obtained for the four-quark 
systems. Finally, in Sect. |3we resume our most important conclusions. 

II. GENERAL FORMULATION OF THE PROBLEM 

The system of two quarks and two antiquarks with the same fiavor can be regarded as 
a system of four identical particles. Each particle carries a SU{2) spin label and a SU{3) 
color label. Both quarks and antiquarks are spin | particles, but whereas a quark color state 
belongs to the SU{3) fundamental representation [3], an antiquark color state is a member of 
the fundamental representation [3]. The four-body wave function is a sum of outer products 



of color, spin and configuration terms 

10) = I Color) I Spin) I i?) (1) 

coupled to yield an antisymmetric wave function with a set of quantum numbers that reflects 
the symmetries of the system. These are the total angular momentum quantum number J, its 
projection J^, and the SU{3) color state G (to avoid confusion with the charge conjugation 
quantum number the color degree of freedom will be labeled as G from now on), which by 
assumption must belong for physical states to the SU{3) color singlet representation. Since 
QCD preserves parity, parity is also a good quantum number. Another relevant quantum 
number to the system under consideration it is the C— parity, C, i.e., the symmetry under 
interchange of quarks and antiquarks. 

To obtain a solution of the four-body Schrodinger equation we eliminate the center of 
mass and use the relative, Jacobi, coordinates r/i, r/2, . . . , tja-i- Then we expand the spatial 
part of the wave-function using the HH basis. In this formalism the Jacobi coordinates 
are replaced by one radial coordinate, the hyperradius p, and a set of {3A — 4) angular 
coordinates Qa- The HH basis functions are eigenfunctions of the hyper spherical part of the 
Laplace operator. An antisymmetric A-body basis functions with total angular momentum 
Ja, Ja: color Ga and C-parity C, are given by. 
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are the hyperradial basis functions, taken to be Laguerre functions. 

{QA\KALAMATAYA-,aA) ^ 3^1gL,M,r,r._,«.(^A) (4) 

are HH functions with hyperspherical angular momentum K = Ka, and orbital angular 
momentum quantum numbers {La, Ma) that belong to well-defined irreducible representa- 
tions (irreps) Fi G r2 . . . G F^ of the permutation group-subgroup chain Si C 1S2 . . . (Z Sa, 
denoted by the Yamanouchi symbol [Fa,!^-!] = [F4, F4- 1, . . . ,Fi]. The dimension of the 
irrep F^ is denoted by |Fm| and A^^Xa-i i^ ^ phase factor [25|. Similarly, the functions 

{sl..s\,gi..gA\SASlGACfA,YA-M ^ X^sUgaTaXa-.^a^'^^-'^^'^'-^^^ ^^^ 

are the symmetrized color-spin basis functions, given in terms of the spin projections (s^) 
and color states (gi) of the particles. The quantum numbers aA, and Pa are used to remove 
the degeneracy of the HH and color-spin states, respectively. For the construction of the 
symmetrized HH basis we shall use the algorithm of Barnea and Novoselsky jSj], which 
utilizes the group of kinematic rotations. For the color-spin subspace, we will present in 
the following subsections a method to transform the standard basis into a symmetrized 
color-spin basis with well defined color and C— parity. The calculation of the Hamiltonian 
matrix-elements between the antisymmetric basis functions, Eq. (j21), is practically the same 
as in the nuclear physics case, replacing isospin by color, and the reader is referred to Barnea 
et. al. W 



A. The construction of symmetrized color-spin states 

The starting point for the construction of color-spin states with well-defined permuta- 
tional symmetry, color and spin quantum numbers, are the single particle states. The spin 
state of the z'th particle, quark or antiquark, is given by |s = | s^ = ±|). Using the Clebsh- 
Gordan coefficients one can successively couple the spin states to construct an A-hodj state 
with well defined total spin 

\S^S|SA-^SA-2■■■S^) = [[[11)11^)]^^ • • ||)]^-^"- • (6) 

Theoretically one can adopt the same procedure for the color states replacing the SU(2) 
Clebsh-Gordan coefficients with the appropriate SU{3) ones. In practice, however, the 
SU{3) coefficients are more involved and we shall therefore follow a different path while 
dealing with color. Each quark (antiquark) color state is one of the three states of the 
triplet representation [3] ([3]) of SU{3). These states are eigenstates of the operators (hi, h2) 
which form the Cartan subalgebra of SU{3). Let us denote by \gi) the color state of the z'th 
particle, \gi) = |1), |2), |3) for quark states and |i), |2), |3) for antiquark states, 

hil^i) = gllgi) ; hsl^i) = g^lgi) . (7) 

The eigenvalues {g},gf) (weights) uniquely label the quark (antiquark) states, in particular 
we can write 

K-i-D) 

l(0,+i)). (8) 

An A-hodj color state is a product of the single particle color states. This state does not 
belong to a well-defined irreducible representation of the color SU{3) group, but instead it 
is characterized by the "color projections" G\ = gl + gl + . . . g\ and G\ = gf + g2 + . ■ ■ g% 
very much the same way as the "m-scheme" states used in nuclear structure calculations 
which have J^ as a good quantum number but not J. From now on we shall use the notation 
Ga for the pair of eigenvalues {G\,G\). The combined color-spin single particle state is 
\si = ^s^gi), and the n-particle states are given by 

\xGs) = m9i)\i92)r^\y3)r'...\i9n)r-'^. (9) 

In the following we shall present a method to transform these states into basis states with 
well-defined permutational symmetry, color and C— parity. This method is based on an 
algorithm by Novoselsky, Katriel and Gilmore (NKG) [25| who devised a recursive method 
to split an invariant vector space into states with well-defined permutational symmetry. 
By now, this algorithm has been used to construct symmetrized states in broad range of 
problems, and therefore we shall not give a detailed account of the method but rather outline 
the main ideas. For a clear description of the method the interested reader is referred to 

The main tool of the NKG algorithm is the transposition class-sum operator. For the 
n-particle permutation group iS„, this operator is the sum over all possible transpositions, 
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It turns out that a simultaneous eigenstate of the mutually commuting transposition class- 
sum operator set {C2[iSn],C2[iS„_i], . . . ,C2[»S2]} has definite symmetry with respect to the 
group-subgroup chain Sn D Sn-i D • • • D ^2 D Si. Accordingly, NKG formulated a 
recursive method that starts with one-particle states and constructs symmetrized A-particle 
states through the following procedure: 

a. Add another particle (let say a total of n). 

b. Identify the subspaces invariant under the action of the permutation group 5„. Create 
the appropriate basis. 

c. In each subspace, evaluate the matrix elements of the transposition class-sum operator 
C2[5„], and diagonalize it. 

d. The Sn symmetry of the resulting states is unambiguously identified through the eigen- 
values of C2[5„]. The eigenvectors of the C2[i5„] matrix are the transformation coeffi- 
cients from the old basis (c.) to the new symmetrized one, these are the coefficients 
of fractional parentage. 



e. Ensure phase consistency (see reference j27[ for a detailed account of this delicate 
point). 

f. li n < A repeat the procedure. 

When this procedure is done, the original states are transformed step by step into the desired 
symmetrized A-particle states. 

In the case under consideration, the n-particle states are labeled by a set of good quantum 
numbers consisting of the total spin Sn and its projection S^ (which we shall suppress), the 
color projection Gn = (G^, G^) and the set of Young diagrams r„, r„_i, . . . , Fs, r2, Fi which 
is equivalent to the Yamanouchi symbol Yn |28|. Using these labels one can construct a 
complete set of states labeled by \SnGnYnl3n), where /5„ is an additional label that removes 
the remaining degeneracies. For single particle, there is no degeneracy and one sees that 

|5iGiFi) = |i^i), (11) 

where Si = ^ and Gi = gi or equivalently {G\,Gl) = {gl,gl). To obtain a two-particle 
state labeled by the sequence of Young diagrams F2, Fi, we first couple the spin states, and 
introduce the notation 

\{SiGiTi;g2)S2G2) = [\SiGiTi)\'^g2)f' (12) 

to denote 2-particle states with total spin 5*2 and total color projection G2 = Gi + g2- Here, 
and in what follows, we have used the short hand notation G2 = Gi + g2 for (^3,6^2) = 
{G\ + glyGf + (72)- The permutation operators transform these states within an invariant 
subspace, each of which is characterized by the total spin 5*2 and the total color projection 
G2 = {Gl, 6*2 )• The linear combinations which belong to well-defined irreps of the symmetry 
group can be written in the form, 

|^2G'2F2Fi/?2) = Yl [SiGiTig2\}S2G2T2P2]\iSiGiTi;g2)S2G2) . (13) 

G\g2;G2=Gi+g2 



Assume that we have constructed a (n — 1)— particle states with well-defined per mutational 
symmetry. Following the 2— particle example we first form n-particle states with well-defined 
total spin and total color projection, 

\[Sn-lGn~lTn-lPn-l', gn)SnGn) = [\Sn-lGn-l^n-lPn-l)\29n)] "• (14) 

Note that G„ = G„-i + dn- The desired, symmetrized, n-particle states can now be written 
in the form 

li^n^n^ nPn) / ^ 

[Sn-lGn-lYn-l/3n-ign\}SnGnYn[3n]\{Sn~lGn-lTn-l(^n-l] 9n)SnGn) ■ (15) 

The transformation coefficients [S'„„iG„„iF„_i/5„_i(7„|}S'jiG'„F„/9„] are the color-spin coef- 
ficients of fractional parentage, cscfps. These coefficients are the lines of an orthogonal 
transformation matrix. Therefore they fulfill the orthogonality and completeness relations. 
The evaluation of the matrix-elements of the transposition class-sum operator, C2[i5n]; for 
the non-symmetrized n-particle states, Eq. ()14|) . can be drastically simplified if we recast it 
in the form 

n-l 
C2[Sn]=C2[Sn^l]+Y,i',n). (16) 

i=l 

The first term on the right-hand side of Eq. (fT^ is diagonal in r„_i. Its explicit value can 
be calculated exploiting the following expression for the eigenvalues of C2[(Sn] |29l |. 

^2[r„] = ^^r,(r,-2z + l), (17) 

i 

where Vi is the number of boxes in the i'th row of r„ and the sum extends over all the rows. 
The matrix elements of the second term between states with well-defined iS„_i symmetry 
can be rewritten as j25i] . 



i=l 

= '5r„_„r;_,5y„_,y._^^^ J2 (r„_iK-2|(n-i,n)|r„_iy::_2)- (18) 

The (n — l,n) matrix elements can be evaluated by expanding the (n — 1)— particle states 
in terms of the cscfps calculated previously and by recoupling the spin states in different 
order, 

\('S'n-iG'„_ir„_i/5„_i; gn)SnGn\\n — 1, n)|(5'„_]^G„_]^r„_i/5„_i; g^SnGn) 

= 2_^ [Sn-2Gn-2Yn-2l3n-2gn-l\]Sn~lGn-lYn~lPn~l\ 

Sn-2Gn-2t3n-2gn-19'„_i 

[Sn-2Gn^2Yn-2Pn-2gn-l\}Sn-lGn~lYn-lPn-l] 
^ V^ J 2 2 ^n,n-l I ) 2 2 '-'","-1 

^^-^ 1 Jn~2 ^n '-'n-l J 1 Jn-2 J-n '-'n-l 
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n.n—l 



{Sn,n-ign-ign\{.ri - l,n)\Sn,n~ig'n-ig'n) ■ (19) 

Since the transposition (n — 1, n) commutes with the [n — 2)— particle operators the matrix 
element is diagonal in the {n — 2)— particle states. The matrix element on the right-hand 
side of Eq. flTUj) is the matrix element between two-particle states after we have decoupled 
them from the rest of the system. This matrix element can easily be evaluated to yield, 

{Sn,n-ign-ign\{n - 1, n)\Sn,n-ign-ign) = -(-l)^"'""''^5n-i,9;^<_i,5n • (20) 

B. SU{3) singlet states and C— parity 

In the previous subsection we outlined a method to construct color-spin states with well 
defined permutational symmetry Ya, spin quantum numbers (SaS^), and color projection 
Ga- Still, the SU{3) color symmetry and C— parity of these states is not established, since 
a state with the same color projection can belong to numerous irreps of SU{3). For example 
the color projection (0, 0) can belong to the singlet representation [0] as well as to the octet 
representation, [8]. The SU{3) symmetry of the A-particle states can be established using 
the quadratic Casimir operator, C2[SU{3)]. Any SU{3) state is an eigenvector of C2[SU{3)] 
with eigenvalue that depends on the specific representation of the state. In particular the 
states with eigenvalue are singlet SU{3) states. The A-particle SU{3) generators 

A 

A^ = E^-' (21) 

1=1 

are symmetric with respect to particle permutations. Consequently the quadratic Casimir, 

8 

C2[SU{3)] = \X^-X^ = \J2KK, (22) 

a=l 

commutes with the permutation group Sa, and the states (IS^iGyiYA/^A)} with good quantum 
numbers (SAGAyA) form an invariant subspace for C2[SU{3)]. The calculation of the matrix- 
elements of the quadratic Casimir can be simplified if we rewrite it as a sum of one-particle 
and two-particle terms 

A A 

C2 [5t/(3)] = i 5^ A,^- X^ + i ^ A^ A,^ . (23) 



The first sum on the right-hand side of Eq. (jzajl is just a sum of the single-particle quadratic 
Casimir operator, and its matrix-elements are just 

A 

{SAGAYAf3A\lYl ^'i ■ ^'i\SAGAYA(3'A) = Sp.,p'^-A , (24) 

since the eigenvalue of the quadratic Casimir for both representations [3] and [3] equals 4/3. 
Due to Schur's Lemma the matrix elements of C2[SU{3)] are independent of the particular 

8 



perinutational symmetry state Ya-i- Exploiting this observation we can recast the second 
sum on the right-hand side of Eq. ()23p into the form 



A 



{SaGaYa/3a\\Y1 ^1 ■ >^'j\SaGaYaP'a) 
^A{A_lll ^ l{SAGArAY:i^,(3A\X'A-i-^A\SAGArAY:i_,(3'A) . (25) 

The Xa-i ■ ^A matrix elements can be evaluated by expanding the A-particle states in terms 
of the cscfps, Eq. (fT3j) . to decouple the last two particles, 

{SaGa^aYa~iPa\Xa-i ■ Xa\SaGa^aYa-iPa) 

SA~lGA~lG\_-ll3A-ll3'^_^gAg'A •S'a-2Ga-2/3a-29a-i9a-i 

y<[SA-lGA-lYA-lPA~igA\}SAGAYAPA][S'A^A^lYA^lPA-l9'A\}SAGAYAPA] 
x[SA-2GA-2YA-2PA-2gA-l\}SA-lGA-lYA-lPA-l] 
x[SA-2GA-2YA~2PA-2gA-l\}SA-lG'A-iYA-lPA-l] 

^{9a9a~i\^a-i ■ ^aWa9'a-i) ■ (26) 



The last term on the right-hand side of Eq. (jzbj) is the sum 

{gAQA-ll^A-l ■ ^aIqWa^i) = '^{9A-l\X'aA-l\9'A-l){9A\X'aA\9'A) , (27) 

a 

which can be easily evaluated using the explicit form of the SU{3) generators for the [3] and 
[3] irreps, Eq. (jH)). 

Turning now to establish states with good C— parity, we note that the C— parity operator, 
C, commutes with Sa, the spin operator and the quadratic Casimir of SU{3). Therefore 
the states {\SaGaYaPa)} with good quantum numbers (SaGaYa) that belong to the singlet 
irrep of SU{3) form an invariant subspace for C. In order to evaluate the matrix-elements 
of C we first express our states in terms of the spin coupled single-particle states, Eq. Q. 
Then we can replace each quark state by the appropriate antiquark state and vice versa. 
The eigenvectors of C with eigenvalues c = ±1 that belong to the subspace of C2[SU{3)] 
eigenstates with eigenvalue '0' are the desired physical color-spin basis functions. 

III. A CONSTITUENT QUARK MODEL 

The charmonium spectra has been thoroughly studied using constituent quark models 
for more than 25 years. As an almost pure coulombic system, models based on coulomb 
plu s confinement interactions have been used to describe most of its spectroscopic properties 
|30| . To address the study of four-quark states we make use of a standard constituent quark 
model inspired in these works and originally applied to the study of the nonstrange baryon 
spectra and the baryon-baryon interaction J3l|. This model has been generalized to all flavor 
sectors giving a reasonable description of the meson spectra [2J], the baryon spectra '5I], and 
the scalar mesons once four-quark configurations were included J20|. Within this model 
hadrons are described as clusters of constituent (massive) quarks. This is based on the 

9 



assumption that the hght-quark constituent quark mass appears because of the spontaneous 
breaking of the original SU{3)l ® SU{3)r chiral symmetry at some momentum scale. In 
this domain of momenta, quarks interact through Goldstone boson exchange potentials. For 
the particular case of heavy quarks, chiral symmetry is explicitly broken and therefore the 
boson exchanges do not contribute to the qq interaction. Explicit expression of the boson 
exchange interacting potentials and a more detailed discussion can be found in Ref . J2^ . 

Beyond the chiral symmetry breaking scale one expects the dynamics being governed 
by QCD perturbative effects. They are taken into account through the one-gluon-exchange 
(OGE) potential J3^. The nonrelativistic reduction of the OGE diagram in QCD for point- 
like quarks presents a contact term that, when not treated perturbatively, leads to collapse 
J33[. This is why one maintains the structure of the OGE, but the 5 function is regularized 
in a suitable way. This regularization is justified based on the finite size of the constituent 
quarks and should be therefore flavor dependent [3^. As a consequence, the OGE reads. 



VocE{n,) = -^a. (A^. ■ A-) |- - ^^^a. ■ a, ^-;:t^\ ^ (28) 

where A*^ are the SU{3) color matrices, a^ is the quark-gluon coupling constant, and ro(/i) = 
rofJ'nn/fJ'ij, where fiij is the reduced mass of the interacting quarks ij {n stands for the light 
quarks u and d) and fo is a parameter to be determined from the data. 

The strong coupling constant, taken to be constant for each flavor sector, has to be scale- 
dependent when describing different flavor sectors 35|. Such an effective scale dependence 



has been related to the typical momentum scale of each flavor sector and assimilated to the 
reduced mass of the system J36|. We use a strong coupling constant given by, 

where Oq, /io and 70 are parameters fitted within a global description of the meson spec- 
tra 0. 

Finally, any model imitating QCD should incorporate confinement. Lattice calculations 
in the quenched approximation derived, for heavy quarks, a confining interaction linearly 
dependent on the interquark distance. The consideration of sea quarks apart from valence 
quarks (unquenched approximation) suggest a screening effect on the potential when increas- 
ing the interquark distance J37|. Creation of hght-quark pairs out of vacuum in between the 
quarks becomes energetically preferable resulting in a complete screening of quark color 
charges at large distances. String breaking has been definitively confirmed through lattice 
calculations J38[ in coincidence with the quite rapid crossover from a linear rising to a flat 
potential well established in SU(2) Yang-Mills theories J39|. A screened potential simulating 
these results can be written as 

VcoN{r,,) = -a, (1 - e-''^'--)(A"'-. ■ A"'-,) . (30) 

At short distances this potential presents a linear behavior with an effective confinement 
strength a = Qc^d^'^i ■ ^'^j), while it becomes constant at large distances. Such screened 
confiningpotentials provide with an explanation to the missing state problem in the baryon 
spectra '401], improve the description of the heavy meson spectra J41|, and justify the devia- 
tion of the meson Regge trajectories from the linear behavior for higher angular momentum 
states 113. 
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IV. RESULTS AND DISCUSSION 

The results for the cccc states have been obtained in the framework of the hyperspherical 
harmonic formahsm explained above up to the maximum value of K within our computa- 
tional capabilities, being K^ax = 20 for positive parity states and K^^ax = 21 for negative 
parity states with the exception of J^*" = l"*"" {K^ax = 18), 2++ and 2+~ {K^ax = 22) and 
2 {Kmax = 23). The two-body problem has been solved using the Numerov algorithm 
(see Ref . \2m for details) and all the model parameters fixed in the description of the qq 
spectra [20|, |2J]. Since the only relevant two-meson decay thresholds for these four-quark 
systems are those formed by two cc mesons, we summarize in Table H] the results obtained 
for the charmonium spectrum in Ref. J2^ compared with the experimental data quoted by 
the Particle Data Group (PDG) |43|. It can be observed that the energy difference in the 
P— wave cc multiplets due to the non-central terms is less than 50 MeV. Although these 
terms are known to play an important role in the description of the light qq states, their 
contribution becomes smaller as the mass of the heavy quark increases J2J] . This allows to 
neglect them on the four-body sector, because they would only provide with a fine tune of 
the final results but making the solution of the four-body problem much more involved and 
time consuming. 

To analyze the stability of these systems against dissociation through strong decay, cccc —>■ 
Mi(cc) + M2(cc), parity (P), C— parity (C), and total angular momentum (J) must be 
preserved. Assuming that the decay takes place in a relative 5— wave for the final state, 
Jiq = Jmi ® Jm2 s-iid P^q = Pmi " Pm2- Siuce the final states are also eigenstates of C— parity 
one has C^g = Cmi ■ C*M2- Once the meson masses have been obtained, the corresponding 
thresholds can be computed simply adding the mesons masses Mi and M2, T(Mi,M2) = 
Ml + M2. In Table [H] we indicate the lowest two- meson decay threshold for each set of 
quantum numbers. Four-quark states will be stable under strong interaction if their total 
energy lies below all possible, and allowed, two-meson thresholds. It is useful to define 

A = M{qiq2q-m) - T{Mi, M^) , (31) 

in such a way that if A > the four-quark system will fall apart into two mesons, while 
A < will indicate that such strong decay is forbidden and therefore the decay, if allowed, 
must be weak or electromagnetic, being its width much narrower. 

Let us first of all analyze the convergence of the expansion in terms of hyperspherical 
harmonics. We show in Fig. ^the variation of the energy of the J^*" = O"*""*" and O""*" cases 
with K, this behavior being similar for the other quantum numbers. From this figure it can 
be seen that the convergence is slow, and the effective potential techniques [23] are unable 
to improve it. In order to obtain a more adequate value for the energy we have extrapolated 
it according to the expression 

E{K) = EiK = oc) + ^^, (32) 

where E{K = 00), a and b are fitted parameters. In Table ITTTl we show the values obtained 
for E{K = cxd) for the two states shown in Fig. ^as a function of the fitting range {Kq, Kf). 
It can be checked that the values obtained for E{K = 00) are stable within ±10 MeV 
for Kq > 10. The values obtained for the b parameter are very similar for all possible 
quantum numbers, being in all cases in the range 0.75 — 0.95, corresponding the lower limit 
to the positive parity states and the upper one to the negative parity states. The origin of the 
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difficulties to obtain a better convergence can be traced back to two different aspects of these 
systems. On the one hand, the HH formahsm is better suited to deal with bound states, 
however, in this particular case we will show that most of the four-quark states are above the 
corresponding two-meson threshold. On the other hand, the color structure of a four-quark 
system is more involved than that of qq or qqq systems J20|, with two different possibilities 
to obtain a color singlet. This makes possible the existence of repulsive contributions due 
to the color operator of the interacting potential, Eqs. (PH|) and (jSHI)- Although the total 
interaction is always confining, these deconfining terms make the convergence much slower. 

We have studied all possible J^*" quantum numbers with L = 0. We show in Table 
|IV| the results obtained for K = Kmax, the maximum value of K calculated, the results 
obtained using the extrapolation of Eq. (jH^ . the corresponding threshold for each set of 
quantum numbers, and the value of A, Eq. (jHT|) . A ffist glance to the results shows that 
there are two sets of quantum numbers, J^^ = 0^^ and J^*" = 2"* , where the four- 
quark configuration is clearly below the corresponding two-meson threshold for dissociation. 
The J^'-' = 1~^~ and J^*" = 2"*"+ states are very close to the threshold and all the other 
quantum numbers are far above the corresponding threshold. Being all possible strong 
decays forbidden they should be narrow states with typical widths of the order of a few 
MeV. It is also interesting to note that the quantum numbers 0^* and 2+^ correspond to 
exotic states, those whose quantum numbers cannot be obtained from a qq configuration, 
and therefore, if experimentally observed, they would be easily distinguished as a clear signal 
of pure non— gg states. 

The energies obtained for the negative parity states are, as expected, much higher than 
those of positive parity and therefore all of them are far above the corresponding thresholds, 
A( J^*^ = J~'") > 500 MeV. Thus, they will fall apart immediately and therefore they should 
be very broad and difficult to detect. Let us also note that the experimental observation 
stating that a unit of angular momentum (parity change) in the qqq and qq sector costs 
approximately 400—500 MeV of excitation energy in all flavor sectors (with some remarkable 
exceptions as the light scalar mesons [20|) has an equivalent in the four-quark sector, being 
in this case the energy difference due to the parity change E{J~^) — E{J^^) k, 800 — 900 
MeV. 

The existence of non— gg signals in the meson spectra has been the subject of an intensive 
debate [lO| . To analyze whether the existence of bound states with exotic quantum numbers 
could be a characteristic feature of the heavy quark sector or it is also present in the light 
sector we show in Fig. |21the value of A as a function of the quark mass for all exotic quantum 
numbers and in Fig. El for the remaining positive parity states. Since each quantum number 
has a different threshold these figures should be interpreted carefully. One should notice that 
the value A = in both figures corresponds to M(gig2g3g4) = T(Mi,M2), and therefore, 
the fact that one state is below the others in these figures do not imply that its total mass 
would be smaller. Let us also note that since the heavy quarks are isoscalar states, the 
flavor wave function of the four heavy-quark states will be completely symmetric with total 
isospin equal to zero. Therefore, one should compare the results obtained in the light-quark 
case with a completely symmetric flavor wave function, i.e., the isotensor states. 

In Fig. Elwe observe how one of the four-quark non-exotic positive parity states, the 2"*"+, 
becomes more bound when the quark mass is decreased, A k, —80 MeV for mq = 313 MeV. 
The l"*"" and 0+^ states, that were slightly above the threshold for rriq = 1752 MeV, increase 
their attraction when the quark mass is increased and only for masses above 3 GeV, close to 
the bottom quark mass, may be bound. With respect to the exotic quantum numbers, we 
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observe in Fig. |21that the and 1 '" are not bound for any value of the quark mass. The 
2"*"" state decreases its binding when the quark mass diminishes and it becomes unbound 
for masses of the order of 500 MeV, the strange quark mass. Only the O"*"" four-quark state 
becomes more deeply bound when the constituent quark mass decreases, and therefore only 
one bound state with exotic quantum numbers would remain in the light-quark sector. 

Taking the experimental mass for the threshold T(Mi,M2) in the light-quark case, one 
can estimate the energy region where these states can be found, being M{0~^^) ^ 900 MeV, 
M(l+-) Ri 1200 MeV, M(l++) Ri 1900 MeV, and M(2++) ^ 1500 MeV. Note that the 
boson exchange potentials, that as explained in Sect. II are present in the description of 
the light-quark sector, have not been considered in this analysis. Although these terms 
do play a role in the description of the four-quark state total energy, it was shown in 
Ref. [1£] that their suppression will not alter the relative order of the four-quark states. 
The four-quark states with quantum numbers O"^"*" and 2"*""*" have been analyzed with the 
complete constituent quark model including boson exchanges using a variational approach 
[20|, obtaining 1004 MeV and 1500 MeV respectively, both in good agreement with our 
calculation and therefore giving confidence on our approach neglecting the boson exchange 
interactions in this particular analysis. The states with exotic quantum numbers would 
be higher in energy, being M{0—) ^ M(l-+) ^ 2900 MeV, M{0+-) ^ 1800 MeV and 
M(2+-) ^ 2100 MeV. 

There are experimental evidences for three states with exotic quantum numbers in the 
light-quark sector. Two of them are isovectors with quantum numbers J^*" = 1 '" named 
7ri(1400) and 7ri(1600), and one isotensor J^^ = 2++, the X(1600) 0. Based on the 
coupling of the 7ri(1400) to the rjir decay channel a possible four-quark structure has been 
suggested [44] while the strong ri'ir coupling of the 7ri(1600) makes it a good candidate for a 
pure hybrid state J45|. The large mass obtained in our analysis for these quantum numbers 
makes doubtful the identification of any of them, and in particular of the 7ri(1400), with a 
pure four-quark state, although a complete calculation with the proper flavor wave function 
is needed before drawing a definitive conclusion |46|. An alternative explanation of the 
7ri(1400) in terms of low-energy rescattering effects has been proposed in Ref. f47] while 
recent reanalysis of experimental data from E852 Collaboration with improved statistics 
show no evidence of any exotic meson with a mass close to 1.6 GeV |48|. Concerning the 
X(1600), being its experimental mass 1600±100 MeV, a possible tetraquark configuration 
seems likely. 

V. SUMMARY 

In this work we have presented for the first time a generalization of the HH formalism 
to study systems made of quarks and antiquarks of the same flavor. For this purpose we 
made use of standard techniques widely used in nuclear physics that have been adapted to 
include the color degree of freedom and to treat explicitly the C— parity. This formalism 
opens the door to an exact study of multiquark systems up to now described by means of 
different techniques, being the variational methods the most standard ones. The particular 
color structure of four-quark systems makes the convergence slow but attainable. 

The formalism has been applied to study the existence of L = four-charm quark states. 
For this purpose we made use of a well-established constituent quark model that properly 
accounts for the charmed meson and baryon spectra. Our results suggest the possible exis- 
tence of three four-quark bound states with quantum numbers 0+^, 2+" and 2++ and masses 
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of the order of 6515, 6648, and 6216 MeV. The two states with exotic quantum numbers, 
clearly below their corresponding two-meson threshold, should present narrow widths and, 
if produced, may be easily detected. 

We have analyzed the variation of our results with the constituent quark mass. In the 
light-quark case only the 0^^ and 2++ quantum numbers remain bound, being the O"*"" the 
lowest one with an energy close to 1800 MeV. The four-quark state with quantum numbers 
1 ^ lies around 2900 MeV, far from the experimental states 7ri(1400) and 7ri(1600), therefore 
supporting a hybrid configuration or a more complicated structure in terms of non-resonant 
structures. A possible description of the X(1600) as a four-quark has also been justified. 

The program we have started for and exact study of multiquark systems by means of 
the HH formalism will be accomplished by implementing the possibility of treating quarks 
of different masses. When this is done we will have at our disposal a powerful method, 
imported from the nuclear physics, to study in an exact way systems made of any number 
of quarks and antiquarks coupled to a color singlet. 
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TABLE I: Charinoniuin spectrum in MeV. Experimental data are taken from PDG |43| . 



(nL) JP^ 


State 


CQM 


Exp. 


(15)0-+ 


Vc{lS) 


2990 


2979.6±1.2 


(15)1- 


J/V^(15) 


3097 


3096.916±0.011 


(1P)0++ 


Xco(l^) 


3443 


3415.19±0.34 


(IP) 1++ 


XclilP) 


3496 


3510.59±0.10 


(IP) 2++ 


Xc2{lP) 


3525 


3556.26±0.11 


(1P)1+- 


hc{lP) 


3507 


3526.21±0.25 


(25)0-+ 


Vc{2S) 


3627 


3654±10 


(25)1- 


V^(25) 


3685 


3686.093±0.034 


{ID)1- 


V'(3770) 


3776 


3770±2.4 


{ID)2- 


^^(3836) 


3790 


3836±13 


(35)1- 


7/^(4040) 


4050 


4040±10 


{2D) 1- 


^^(4160) 


4104 


4159±20 



TABLE IL Lowest 5— wave two-meson thresholds (MeV) for all J quantum numbers. 



JP^ 


Ml Ma 


T(Mi,M2) 


0++ 


7?e(15) r,,(15) 


5980 


0+- 


XciilP) hcilP) 


7003 


1++ 


J/V'(15) J/^(15) 


6194 


1+- 


J/V^(15) 7?e(15) 


6087 


2++ 


J/V(15) J/V'(15) 


6194 


2+- 


rjcilS) ^^(3836) 


6780 


0-+ 


r/,(15) XcoilP) 


6433 


0— 


J^IS) XciilP) 


6593 


1-+ 


Vc{lS) Xci(lP) 


6486 


1— 


r?,(15) hcilP) 


6497 


2-+ 


VcilS) Xc2{lP) 


6515 


2— 


J/V(15) XciilP) 


6593 
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TABLE III: E{K 
JP^ = 0-+. 



oo) (MeV) as a function of the fitting range {KQ,Kf) for J^*-^ = 0^+ and 



{Ko,Kf) 


^(is: = oo)[o++] 


{Ko,Kf) 


EiK = oo)[0-+] 


(2,20) 


5831 


(3,21) 


7018 


(4,20) 


5947 


(5,21) 


7007 


(6,20) 


5997 


(7,21) 


7004 


(8,20) 


6017 


(9,21) 


7003 


(10,20) 


6030 


(11,21) 


7000 


(12,20) 


6038 


(13,21) 


7000 


(14,20) 


6041 


(15,21) 


6993 


(16,20) 


6047 


(17,21) 


6990 



TABLE IV: cccc masses obtained for the niaximum value of K computed, E{Kmax)^ and using the 
extrapolation of Eq. (|32|) . E{K = cxd), compared with the corresponding threshold for each set of 
quantum numbers, T(Mi,M2), as given in Table ITTI The value of A for each state is also given. 
All energies are in MeV. 



JP^ 


^yJ^max 


) E{K = oo) 


T{Mi,M2) 


A 


0++ 


6115 


6038 


5980 


+58 


0+- 


6606 


6515 


7003 


-488 


1++ 


6609 


6530 


6194 


+336 


1+- 


6176 


6101 


6087 


+14 


2++ 


6216 


6172 


6194 


-22 


2+- 


6648 


6586 


6780 


-194 


0-+ 


7051 


6993 


6433 


+560 


0— 


7362 


7276 


6593 


+683 


1-+ 


7363 


7275 


6486 


+789 


1-- 


7052 


6998 


6497 


+501 


2-+ 


7055 


7002 


6515 


+487 


2— 


7357 


7278 


6593 


+685 
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FIG. 1: Energy of the 0^+ (circles) and ^ (crosses) four-quark states as a function of K. 
The sohd hne corresponds to the extrapolation used for the 0^"'' case while the long-dashed one 
corresponds to the '" state. The thresholds for each state are denoted by a dashed line (0^^) 
and a dotted-dashed line (0 ^). 
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FIG. 2: Variation of the energy of the four-quark states with exotic quantum numbers with the 
quark mass. The sohd line corresponds to the quantum numbers J^^ = 0^ , dashed to 2^ , 
long-dashed to , and dashed-dotted to 1 ^. The vertical lines correspond to the light and 
charm quark masses, respectively. 
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FIG. 3: Variation of the energy of the non-exotic positive parity states with the quark mass. The 
grey box corresponds to the ±10 MeV uncertainty in the extrapolation. The sohd hne corresponds 
to the quantum numbers J^*-" = 2++, dashed to 1"^ , long-dashed to 0^+, and dashed-dotted to 
1++. The vertical lines correspond to the light and charm quark masses, respectively. 
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